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1/^ ' Abstract 

■ We provide a simple mathematical description of the exchange of en- 

ergy between two fluids in an expanding Friedmann universe with zero 
^ , spatial curvature. The evolution can be reduced to a single non-linear dif- 

rC) ' ferential equation which we solve in physically relevant cases and provide 

\^ I an analysis of all the possible evolutions. Particular power-law solutions 

f^ . exist for the expansion scale factor and are attractors at late times un- 

^^ ' der particular conditions. We show how a number of problems studied 

^^ , in the literature, such as cosmological vacuum energy decay, particle an- 

nihilation, and the evolution of a population of evaporating black holes, 
correspond to simple particular cases of our model. In all cases we can 
CJ> , determine the effects of the energy transfer on the expansion scale factor. 

O^' We also consider the situation in the presence of 'anti-decaying' fluids and 

t I ' so called 'phantom' fluids which violate the dominant energy conditions. 

1 Introduction 

X 

$H ' There are many cosmological situations where the transfer of energy between 

two fluids is important. The interaction between matter and radiation ^, the 
decay of massive particles into radiation i^, matter creation |3], the forma- 
tion and evaporation of primordial black holes |3], the annihilations of particle- 
antiparticle pairs 0, particle or string production [HlEl, inflaton decay |H| and 
the decay of some scalar field or vacuum energy [TO], are all particular ex- 
amples which have been studied in general-relativistic cosmology. The situation 
in Brans-Dicke cosmology has also been investigated QT], as have the cases of 
two arbitrary interacting fluids ■12" and more than two interacting fluids |13) . 
In some cases, as in the example of accreting and evaporating black holes, there 
will be a two-way transfer of energy occurring as, say, a spectrum of radiation 
inhomogeneities collapse under their self-gravity in the early universe to form 
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a population of primordial black holes but the products of the Hawking evap- 
oration of the black holes adds to the cosmological population of interacting 
relativistic particles 0] . The different studies of these particular situations have 
often identified the existence of special power-law scaling solutions. In this pa- 
per we consider a general problem of this sort, describe its general behaviour, 
relate it to the existence of special power-law solutions, and describe its general 
solution succinctly in terms of the parameters defining the energy exchanges. 
The examples in the literature can then be shown to be particular examples of 
these solutions and the conditions for their stability are made clear. 

We will consider the mutual exchange of energy between two fluids at rates 
that are proportional to a linear combination of their individual densities and 
the expansion rate of the universe. In the absence of any interaction the fluids 
reduce to two separate perfect fluids. 

2 Decaying Fluids 

Consider a flat Fricdmann Robertson Walker (FRW) universe with expansion 
scale factor a(t) containing two fluids with equations of state 

P = (7-l)p, 
Pi = (r-l)pi, 

where the 7 and T are constants, and the evolution of the Hubble parameter 
H = aj a is governed by the Friedmann equation 

31/2 ^p + pi, (1) 

where SttG = 1. Assume that the two fluids exchange energy but the total 
energy is conserved so that 

pi-h3iirpi = -fiHpi + apH, (2) 

p + 3i/7p = [iHpi - apH, (3) 

where a and /3 are constants parametrising the energy exchanges between the 
two fluids. Generalisations of this simple cosmology to spatially curved or 
anisotropic universes can be made in an obvious way if required [3 ^]. In 
an expanding universe {H > 0) this scenario corresponds to p and pi 'decaying' 
into each other in proportion to their energy densities if a and /? are positive. 
The degenerate case 7 = F can be seen to be trivially equivalent to the standard 
scenario without energy exchange, by considering the fluid P2 = P + Pi- 

Using the last three equations we can eliminate the densities to obtain a 
single master equation for the Hubble expansion, H{t) : 

H + HH{a + /3 -F 37 + 3r) + ^H'^{aT + Pj + 3F7) = 0. (4) 

Let us rewrite eq. (0)) as 

ii + AHH + BH^ = (5) 



with 

A = a + /? + 37 + 3r, 
B = haV + p-f + iT-i). 

This equation is a special case of the more general differential equation^ consid- 
ered by Chimento |15j . In |15j Chimento investigates the mathematical structure 
of this equation by showing that it has a form invariance, which is subsequently 
used to find solutions. 

Simple self-similar solutions to eq. lO exist with 

where, for h ^ Q, 

2- Ah + Bh^ = 0, 

and there are two non-trivial solutions, H+{t) and H-{t), with h values 



These real power-law solutions for H{t) exist iff A^ ^ SB. For a,/3,7,r ^ 
and 7 7^ r, this inequality is always satisfied. We can see this by defining 

B_ ^ 3(ar + /?7 + 3r7) 
~ A^ 2(a + ^ + 37 + 3r)2 ' ^ > 

so that A'^ ^ 85 iff 5 ^ 1/8. We see that the denominator in equation IJHJ) 
is always positive, so 5 is always non-singular and positive for finite and semi- 
definite positive values of a, /3, 7 and F. It can also be seen that (5 — > as either 
a or /3 ^ 00 . The maximum value of 5 must therefore occur at finite values of 
a and (3. If this maximum exists when a and (3 are both non-zero then there 
must exist a point at which 

Using © we can see that this condition is never met, so the maximum value of 
6 must exist when a = 0, when /3 = 0, or when a = f3 = Q. For a = Q we will 
have a maximum at non-zero (3 when 

which occurs iff /? = 3(7 — F), for a, /?, 7, F ^ and 7 7^ F. Similarly, the 
maximum can occur at non-zero a when /? = iff a = 3(F — 7). We can choose, 

^y + afy + f3 J fdy + jf = 0, where y = y{x), f = f{y) and ovcrdots denote differentiation 
with respect to x. a, /B and 7 are constants. 



without loss of generality, F > 7 so that the maximum value of 5 occurs when 
a = 3(r — 7) and /3 = 0, and wc have the conclusion that 

O 

for all 7, r ;s and 7 7^ r. 

Having established the existence of the power-law solutions ((TJ , we can now 
show that they behave as attractors of the general solution by solving Q . For 
A^ > 8B we find the solution 

where ci and C2 are constants. This solution to ij^l was previously found by 
Chimento in 1151. As a ^ 00, we then have 



^2 _^ ^-{A-VA^-8B)/2 

and, as a ^ 0, 

J/-2 _^ ^-{A+VA^^8B)/2 

These two equations can be integrated to obtain 

^^^^iA±VA^^8B)/2B^ (10) 

which are the power-law solutions 10, found earlier in eq. ©. By integrating 
(O we can show explicitly the existence of the above power-law attractors, and 
the smooth evolution of a between them. It is possible to integrate © to get 
a solution in terms of t and the hypergeometric function 2Fi{a,b;c;x). An 
expression in terms of more transp arent f unctions can be found by defining a 
new time coordinate dr = a-l^+v-^ -8B)/4^^ ^^^ integrating in terms of r. 
This gives the solution 



a (X e 



V^ir-r„)^^ _ gV^^(A2-8B)(r-ro)^-2/^(A2-8S)^ (•JJ^) 



This is the same form for the evolution of a that was found by Chimento and 
Lazkoz in their investigation of phantom fluids in k-essence jl6| . It can be seen 
from this expression that a — > as a ~ gv^v (-^ -8B)(r-ro) .^^jj^gj-^ j- _^ —00. 
In terms of the coordinate t, this corresponds to the solution a_ above. As 
a — > 00, the solution smoothly approaches a ^ (t — tq)^^^^^^~^^ as t — > tq, 
which corresponds to the solution a+. 

We have now shown that the two power-law solutions 1)10(1 exist and for all 
a, /3, 7, F > and are the attractors of the smoothly evolving general solution 
at late and early times when A^ > 8B. 

It remains to investigate the limiting case A"^ = 8B. The exact solution to 
equation Q when A"^ = 8B is jTS] 

H^ = a-^^^{c3 + Ci\na), (12) 



where C3 and C4 are constants. For C4 — this solution corresponds to power- 
law expansion described by the degenerate case where a^ = a-. For C4 7^ this 
solution is more complicated and is bounded by a = g-c3/c4 -ftrhjigt approaching 
H^ ^ a~^/^ In a as a ^ or 00, which does not describe a power-law behaviour. 
An illustrative special exact solution to equation Q) exists when B = A^ /9, 
as was shown by Chimento JS]. In this case @ can be linearised to ^/i = by 
the substitution H = {3/A)tJ;/il;. Hence, for this special value of B 

^^ 3(c,+ 2ce^ (13) 

where C5, cg are constants of integration. This expression can be integrated to 

A 

— ] ' =l + C5t + C6t^, (14) 

a-oj 

where oq is a constant and the early and late time behaviour is clear and has 
the same form as the power-law solutions (0) when B — A^ /9. 

3 Evolution of the energy densities 

The conservation equations Q and l|2Jl can be used to construct the second- 
order differential equation 

[L+A^ + 2B = (15) 

P P 

where A and B are defined as before and primes denote differentiation with 
respect to the variable 77 = In a. This equation can be solved for p and the 
corresponding solution for pi can then be found from Q. Substituting these 
solutions into the Friedmann equation (QJ gives, for A^ > B, the solution ^ 
that was previously found by solving the master equation Q . 

The advantage of considering the evolution of p directly is that a particularly 
interesting behaviour can be observed in the evolution of the ratio p/ pi for the 
self-similar solutions (jHJl . To find this behaviour we first note that a solution to 
equation (|15|l is given by 

N 

p^ Poa 



where po is a constant and 2N ~ —A ± \/ A^ — 8B. Substituting this into 
equation Q gives the corresponding solution for pi 

Pi = pioa 

where pio — {N + 3"/ + ajp^/ j3 is constant. These solutions for p and pi, 
when substituted into the Friedmann equation |^ , correspond to the self-similar 
solutions for H given by ©. It is immediately apparent that p and pi evolve 
at the same rate and so the ratio p/pi is a constant quantity 

p p 



px (A^ + 37 + a) 



during a period described by the power-law evolution pU|l . It is this constant 
ratio in the energy density of two fluids with different barotropic indices 7 and 
r that has been used by a number of authors in an attempt to alleviate the 
coincidence problem concerning the the present-day values of the vacuum and 
matter energy-densities [51ll0|. 

4 Three Examples 

The exact solutions found in the last sections provide us with extensions of the 
analysis of several cosmological problems that have been studied in the past, 
which can be defined by particular choices of the two parameters A and B. As 
we have seen, the overall dynamical behaviour is determined by the behaviour 
of the combination 6 = B/A^. 

4.1 Particle-antiparticle annihilation 

Consider the problem of the long-term evolution of a universe containing equal 
numbers of electron-positron pairs jS] . If these particles are assumed to be the 
lightest massive charged leptons then they cannot decay, and can only disappear 
by means of the mutual annihilations e~e"'" -^ 27. Page and McKee set up a 
model for the e~e^ annihilation into radiation that corresponds to taking the 
special case a = 0,/3 > 0,r = 1,7 = 4/3 in equations ^ and © and the 
definition of Ppm by Page and McKee is given in terms of our /3 by /3 = 
3/3pAf /(2 — Ppm)- They find the power-law solution with h = h- = 2/(/3-|- 3) = 
(2 — l3pM)li which reduces to the usual dust FRW model when (i — f3pM = 
and there is no annihilation into radiation. The effect of the annihilations is to 
push the expansion away from the dust-dominated form with a — t^'^ towards 
the radiation dominated evolution with a = t^/^. The other power-law solution 
corresponds to the pure radiation case with h — h+ — 1/2. We can verify that 
this power-law solution is an attractor by evaluating S, since for the e~e+ —^ 27 
annihilation /3pa/ = (13 - VT05)/8 = 0.3441 so S = B/A^ = 0.1247 < 1/8. 

4.2 Primordial black-hole evolution 

A more complicated energy exchange problem was formulated by Barrow, Copcland 
and Liddle ^ who consider the problem of a power-law mass spectrum of 
primordial black holes forming in the early universe and then evolving un- 
der the effects of Hawking evaporation of the part of the mass spectrum with 
Hawking lifetimes less than the expansion age. This has two effects. The 
radiation background is supplemented by input from the black-hole evapo- 
ration products and the fall in the total black-hole density goes faster than 
the adiabatic pbh oc a~^ that occurs in the absence of decays because the 
black hole population is a pressureless gas to a very good approximation, since 
p/p ~ w^ ~ T/Mbh ~ {rripi / M}jh){tpi / tY^"^ ~ for masses less than the Planck 
mass nipi at times greater than the Planck time tpi. Accretion of background 
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radiation in the radiation era of the universe by the black holes could be in- 
cluded, but is negligible. This corresponds to our model in the special case 
r = 1, 7 = 4/3,a = and 

8 — n 
where the initial number density spectrum of black holes with masses between 
m and m + 5ni at time t is given by N{m, t) ex to^" and n > 2. 

A power-law solution was found in ref. .2j with /i = (8 — n)/9, so long as 
2 < n < 7/2, and the black hole evaporations have a significant effect on the 
expansion rate of the universe during the radiation era. We have 

A = 

B = 

and so 

~ A2 (25-27i)2- 

We see that the allowed range of n S (2, 7/2) corresponds to J £ (6/49, 1/8) and 
the expansion scale factor evolves as a c>c i(*~")/9. The n = 2 limit corresponds 
to a pure dust-dominated expansion, with a ex t^'^, while the n — 7/2 limit 
corresponds to a pure radiation-dominated evolution, with a ex i^'^. Again we 
see that the power-law solution is an attractor for the general solution with n 
in this range. When the expansion of the universe becomes dominated by cold 
dark matter, with pcdm oc a'~^ , a power-law scaling solution no longer exists 
because the radiation products from the black-hole evaporations now make a 
negligible contribution to the total density of the universe, which is dominated 
by Pcdm > Pbh >> P-y, and a ex t'^'^ becomes the attractor for the evolution of 
the expansion scale factor. 

4.3 Vacuum decay 

The cosmological evolution created by the decay of a vacuum stress (pi — py) 
into equilibrium radiation was considered by Freese et. al. and many other 
authors ^HI ■ It is described by a special case of our equations (O and Q with 
r = 0, 7 = 4/3, a = and /3 > 0. It represents the decay of a scalar field stress 
with p w —p into radiation. In this case we have 

A = /3-t-4, 

B = 2/3, 
, _ B 2/3 



^2 (13 + 4)2 



(16) 



with ft,+ = 1/2 and /i_ = 2//3. We see that the first of these corresponds to the 
degenerate situation with pure radiation. The second solution has a ex P'^ and 



requires /3 > 3 if the evolution of the universe is to have a matter-dominated era 
following a radiation era. As the value (3 increases, the dominance of the vacuum 
contribution slows the expansion whereas in the limit j3 ^> Q the expansion rate 
increases without bound and the dynamics approach the usual vacuum-energy 
dominated de Sitter expansion with a oc exp(i.yp^/3). Again we see that this 
simple solution can be generalised by using the full analysis provided above. We 
see from (|16|) that we always have 5 ^ 1/8 with the maximum of 5 achieved 
when /3 = 4. The solution l|13(l - (|14|l arises for 6 = 1/9 which occurs when 
/3 = 2 or /3 = 8. 

5 Anti-Decaying Fluids 

It was shown in section |21 that for a, /3, 7, T ^ and 7 7^ F the maximum 
value that 5 can take is 1/8. If we relax these assumptions, then 6 can take 
values greater than 1/8 and the qualitative character of the solutions to Q) is 
significantly altered. We will begin by investigating the conditions required for 
5 > 1/8. 

In the previous section it was shown that for a, /3, 7, F ^ and 7 7^ F the 
only point at which 6 = 1/8 is at a = 3(F — 7) and (3 = 0, and at all other 
points in this parameter range we have S < 1/8. It can be seen from that 
S = 1/8 when 

a = (V3(F-7)± ^^y ^ 

and that the first derivatives of 6 are non-zero at any point where this condition 
is satisfied. These values of a therefore separate regions where S < 1/8 from 
those where 5 > 1/8. It can also be seen that S > 1/8 only if a > and 
/3 < 0. These conditions correspond to the fluid p decaying and the fluid pi 
anti-decaying. (By 'anti-decaying' we mean gaining energy in proportion to its 
energy density, instead of losing it). An example of an anti-decaying fluid is a 
ghost field which radiates away energy; here the energy density of the ghost is 
negative, so a negative value of /3 is required for the radiation to carry away 
energy. 

For 6 > 1/8 the exact solution to equation © is 

H^=a'^^^cosf-^/8B~A^\iia\, (17) 

where integration constants have been rescaled into a and H. Again, this equa- 
tion is difficult to solve in term of the coordinate t. By introducing the new 
coordinate dr = a'^^^dt, we get 

din a 

dr ™"V2 

which can be integrated to obtain a closed form for the expansion scale factor: 



= cosi -^/8B - A^lna] , 



4 V8B~^ 

a(T) = exp < , = am ; (r — tq] 

^ ' 1 V8B - A2 I 4 ^ "^ 




(18) 




Figure 1: The Jacobi amplitude am(T|2) 

where tq is constant and ani(a|6) is the Jacobi aniphtude, shown in figure ^ 
The form of a{T) in (|18|l is an always-positive oscillatory function of the time r, 
with constant amplitude. The corresponding solution in terms of the coordinate 
t will therefore also be oscillatory with constant amplitude. 

Whilst in the previous section we found that, for S < 1/8, the scale factor 
evolves as a smooth function with early and late-time power-law behaviour, we 
have found for S > 1/8 a substantially different behaviour. The scale factor now 
oscillates in time and does not display the simple power-law behaviour found in 
the 6 ^ 1/8 situations. 

6 Phantom Fluids 

We have so far only discussed the cases where 7, F ^ and 7 7^ F. This 
assumption is useful as it means that 6 is non-singular in the parameter range 
a, P ^ 0, for which it was shown in section |21 that the maximum value of 6 is 
1/8. This result was subsequently used in section |S1 to show that there exists 
a parameter range with a > and /3 < for which S > 1/8. In this section 
we will relax the positive semi-definite assumption on the parameters 7 and 
F, extending the analysis we have so far performed to the case of so called 
'phantom' fluids. We begin by showing that in the parameter range a, /? > 
there exist no points at which S — > +00. For this to occur we would require the 
simultaneous satisfaction of the conditions 



aT + (3-f + 3F7 > 



and 



a -I- /3 + 37 + 3F = 0. 



Using the second of these conditions, we can cUminate a in the first to find 
/3(7 - r) - 3r2 > or (3<~- — - «; 0, 

(r-7) 

where we still assume F > 7, without loss of generality. Similarly, we can obtain 
for a the expression 

These two inequalities show that 6 — > +00 can only occur in the parameter 
space a > and /3 < 0. Therefore, in the range a, (3 > the only singularities 
in i5 that can occur are those in which (5 ^ —00. In this case we can again show, 
using the arguments in section [3 that the maximum value of S when a, /? ^ 
is 1/8. The argument showing the existence of a region where S > 1/8 in the 
range a > and /3 < now follows in exactly the same way as for the 7, F > 
case, given in sectional 

The form of the solutions in the regions where S < 1/8 and S > 1/8 are the 
same as in the non-phantom case, and are given by equations l|njl and (|17|l . 



7 Discussion 

We have determined the general solution of a simple model with the exchange 
of energy between two fluids in an expanding Friedmann universe of zero spatial 
curvature. The total energy of the exchange is conserved and the model allows 
energy inputs and outflows proportional to the densities of the two fluids. A 
number of simple examples of this sort already exist, such as particle decays or 
particle-antiparticle annihilations into radiation, particle production, the evap- 
oration of a population of primordial black holes, the decay of a cosmological 
vacuum or cosmological 'constant', and energy exchanges between quintessence 
and ordinary matter or radiation. However, these examples are restricted to 
one-way energy exchange and do not prove that the scaling solutions that they 
employ are attractors for the general solution. We have established the exis- 
tence and form of simple power-law solutions for the expansion scale factor in 
the case of two-way energy exchange between fluids and determined that they 
are attractors for the late-time evolution in situations that are usually regarded 
as generic. If we allow one fluid to be anti-decaying then we can move into a 
domain where these power-law solutions are no longer attractors. Again, we find 
the general behaviour for these cosmologies. These solutions provide a simple 
model for the study of a wide range of energy exchange problems in cosmology 
and also reveal the conditions under which power-law solutions previously used 
to solve some of these problems are stable attractors. They provide a simple 
model for many future studies of a variety of interacting fluid cosmologies. 
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